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ABSTRACT 

It has become clear in recent years that globular clusters are not simple stellar populations, but 
may host chemically distinct sub-populations, typically with an enhanced helium abundance. 
These helium-rich populations can make up a substantial fraction of all cluster stars. 

One of the proposed formation channels for blue straggler stars is the physical collision 
and merger of two stars. In the context of multiple populations, collisions between stars with 
different helium abundances should occur and contribute to the observed blue straggler popu- 
lation. This will affect the predicted blue straggler colour and luminosity function. 

We quantify this effect by calculating models of mergers resulting from collisions be- 
tween stars with different helium abundances and using these models to model a merger pop- 
ulation. We then compare these results to four observed clusters, NGC 185 1, NGC 2808, NGC 
5634 and NGC 6093. 

As in previous studies our models deviate from the observations, particularly in the colour 
distributions. However, our results are consistent with observations of multiple populations in 
these clusters. In NGC 2808, our best fitting models include normal and helium enhanced 
populations, in agreement with helium enhancement inferred in this cluster. The other three 
clusters show better agreement with models that do not include helium enhancement. We 
discuss future prospects to improve the modelling of blue straggler populations and the role 
that the models we present here can play in such a study. 

Key words: stars: blue stragglers, evolution, luminosity function - globular clusters: general 



1 INTRODUCTION 

Blue stragglers are stars that appear on the extension of the main 
sequence in both open and globular clusters and have been stud- 
ied extensively since their original discovery by Sandage (1953). 
Their existence cannot be explained by conventional stellar evolu- 
tion models. Two formation channels have been discussed exten- 
sively in the literature, one using stellar collisions (Lombardi et al. 
[T995l|Sills et al.|1997||2005HGlebbeek & Pols|2008) and one using 

binary evolution dChen & Han|2004| 12008 2009). So far observa- 
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tions do not clearly indicate which of these formation mechanisms 
is the dominant one, or even if there is a dominant formation mech- 
anism ( Knigge et al. 2009 1. Neither formation mechanism has been 
able to completely account for the properties of the observed blue 
straggler populations, in particular the colour and magnitude distri- 
butions. 

Star to star abundance variations in globular clusters, origi- 
nally known as "globular cluster abundance anomalies", have been 
known for a long time (see Gratton et al. 2 004| for a recent review). 
The observed abundance patterns indicate that at least some of the 
globular cluster stars (up to 50% of the present day population) con- 
tain material that has undergone partial hydrogen burning at high 
temperatures. The most famous example is the so-called oxygen- 
sodium anti-correlation. To date, abundance variations have been 
observed in all clusters where the observational data has sufficient 



resolution to detect them. Despite work by several groups, the ori- 
gin of these abundance variations is still unclear. The leading mod- 
els propose that material has been processed by AGB stars (Ventura| 
|et al.|2001| |2002| |Ventura & D'Antona|2005) or rapidly rotating 
massive stars (Decressin et al. 2007) before being returned to the 
interstellar medium to form new stars. An alternative scenario in- 
volves non-conservative mass transfer in a binary system ( |de Mink] 
|et al.| [2009). Since the main product of hydrogen burning is he- 
lium, this polluted material is also expected to be relatively helium 
rich and we may therefore expect to find a helium enhanced sub- 
population among cluster stars. 

The first evidence for the existence of multiple populations in 
gobular clusters comes from the discovery of a split giant branch in 
at Cen (Pa ncino et al.|2 000 . It was later shown that the sub giant 
branch JFerraro et al.f 2004 1 and the main sequence JBedin et al.| 
|2004| l are likewise split in multiple sequences, which is interpreted 
as a spread in helium abundance ( |Piotto et al.||2005) . |D'Antona| 
|et al. | (2005 ) predicted the existence of multiple populations in NGC 
2808 based on the observed width of the main sequence band. De- 
tailed observations by Piotto et al. (2007 ) found that the main se- 
quence of NGC 2808 indeed splits into three distinct sequences. 
Milon e et al.| (2008 ) reported the existence of a double sub-giant 
branch in NGC 1851. 

If collisions contribute to blue straggler formation and the dif- 
ferent sub-populations have a comparable number of stars, then col- 
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Table 1. Parameters for the different model sets 



Model set 


Y 


'coll 

Myr 




M,,M 2 
M 


A 


0.24 


6000- 


12000 


0.4, 0.6, 0.8 


B 


0.32 


6000- 


12000 


0.4, 0.6, 0.8 


C 


0.40 


6000- 


12000 


0.4, 0.6, 0.8 


D 


0.24; 0.32 


6000- 


12000 


0.4, 0.6, 0.8 


E 


0.24; 0.32; 0.40 


6000- 


12000 


0.4, 0.6, 0.8 



lisions between stars from a helium-rich and a helium-normal pop- 
ulation must contribute to the observed blue straggler population. 
A blue straggler that is formed from a collision between a helium 
normal and a helium-rich star will have a higher helium abundance 
in its interior than a blue straggler formed by a collision between 
equivalent stars with the same helium abundance. This has the ef- 
fect of making the star brighter, even if the material is not mixed 
efficiently 1 Gleb beek et al.| [2008). Therefore, collisions involving 
stars from different populations should be taken into consideration 
when constructing the blue straggler colour and luminosity func- 
tions in clusters with multiple populations. 

Here we explore the results of such a calculation. 



sions between stars with Y = 0.24 and Y Q = 0.32 and is a superset 
of models A and B. Similarly model set E is a superset of model 
sets A, B and C that also contains the cross-collisions involving 
stars with Y = 0.40. 

Model sets A, B and D have been calculated for Z = 0.0003 
as well as Z = 0.001. 



2.1 MMAS 

The structure of our collision products was calculated using the 
Make Me A Star (MMAS) code by |Lombardi^nLl ( (2002l >. MMAS 
approximates the structure of a collision product using an algorithm 
known as entropy sorting. 

The idea behind this algorithm is that the quantity A = p/p 5/3 , 
which is related to the thermodynamic entropy, increases (nearly) 
monotonically inside stars from centre to surface. In the absence of 
strong shocks, as is the case for low-velocity collisions, A is con- 
served in the fluid elements from both colliding stars. The structure 
of the collision product can then be approximated by sorting the 
mass shells of both parent stars in order of increasing A. As a result 
of stellar evolution A decreases in the core, so that the core of the 
collision product is most likely to resemble the core of the most 
evolved parent star ( |Sills et al.|1997)|Glebbeek & Pols|200"8) . 



2 COMPUTATIONAL SET-UP 



2.2 The STARS code 



We use the STARS code ( §2. 2} to calculate the evolution tracks of 
both the progenitor stars and the collision products. The structure 
of the collision products themselves is calculated using the MMAS 
code (SfTT}. 

The parameter space for our models consists of the masses of 
the two colliding stars, Mi and Mi, the time of collision / co n and the 
helium abundance of the colliding stars, K] and Y 2 . In principle the 
impact parameter for the collision spans another dimension in pa- 
rameter space, but we will limit ourselves to non-rotating collision 
models here. In essence, we assume that the collision product has 
an efficient way to lose excess angular momentum {e.g. by the ac- 
tion of a magnetic field) and that rotational mixing is unimportant. 
With this assumption, the effect of a non-zero impact parameter on 
the structure of the collision product is small (Sills et al. 1997) and 
can be ignored. 

We calculated a grid of models that maps out the relevant 
portion of parameter space for blue straggler formation in glob- 
ular clusters. Our computational grid covers a mass range of 
0.4 . . . 0.8M G for both stars. We will refer to the most massive star 
as the primary and the least massive one as the secondary. To ac- 
count for the evolution of the parent stars as well as the evolution of 
the collision product since time of collision the collision time is var- 
ied between 6 000 Myr and 12 000 Myr in steps of 1000 Myr. The 
heavy element content Z is set to Z = 0.001 for most of our mod- 
els and abundances are scaled to solar values. The initial (ZAMS) 
helium abundance in the stars is set to Y = Yq + 2Z, where Fq 
is 0.24,0.32 or 0.40. This parameter range was chosen to cover 
both normal helium abundances (Y = 0.24) and the most extreme 
helium abundance proposed in the literature (Yq = 0.40, in NGC 
2808, |D'Antona et al.|2005]|Piotto et al.|2007) . The initial hydro- 
gen abundance X = 1 — Y — Z and the primary and secondary are 
allowed to have a different Yq. The range of parameters in our grid 
is summarised in TableQ] 

The model sets A, B and C all involve just collisions between 
stars of the same helium content. Model set D also includes colli- 



Our evolutionary models are calculated using a version of Eggle- 
ton's stellar evolution code (Eggleton 1971 ; Pols et al. 1995 ), here- 
after STARS. The STARS code solves the equations of stellar struc- 
ture and the nuclear energy generation rate simultaneously on an 
adaptive non-Lagrangian non-Eulerian ("Eggletonian") grid ( |Stan-| 
|cliffe|2006j . Since |Glebbeek et al.|f2008) we have update the nu- 
clear reaction rates to the recommended values from ( |Angulo et al.| 
|1999| >, with the exception of the 14 N(p, y) 15 reaction, for which we 
use the recommended rate from |Herwig et al.| {2006} and 



Formi- 



|cola et aL| (12004). We use the opacity tables of Eldridge & Tout 



( 2004 1, which combine the OPAL opacities from Iglesias & Rogers 
(1996), the low temperature molecular opacities from Alexander 



& Fe rgusonM1994) , electron scattering opacities fromjBuchler & 



Yueh| q 1976 ) and the conductive opacities from Iben ( 1975]. The 
assumed heavy-element composition is scaled to the solar mix- 
ture of (Anders & Grevesse|1989|>. Chem ical mixing due to con- 
vection (|Bohm-Vitense||1958| |Eggleton||1972) and thermohaline 
mixing (Kippenhahn et al. 1980; Stancliffe et al. 2007) is taken 
into account. All models are computed with a mixing-length ra- 
tio l/Hp = 2.0. As in |Glebbeek & Pols|p008) , we have neglected 
convective overshooting in all models described here. Its effect on 
the evolution of stars in the mass range of our collision products 
(0.8M o - 1.6M Q ) is negligible (Pols et al.|1998] >. 

On the giant branch and the early AGB we adopt the Reimers- 
like mass loss prescription from Schroder & Cuntz ( 2005 , 2007) , 



M 



n 



LR 
~M 



V4000K/ 



1 + 



4300g 



(1) 



Because we cannot calculate through the core helium flash, 
we stop the code at the onset of the flash and construct a "zero-age 
horizontal branch" model with the correct total mass, core mass and 
composition. This construction is performed by pseudo-evolving a 
low mass core burning helium star. Composition changes due to 
helium burning are disabled during this process, but hydrogen is 
allowed to burn normally so that the helium core can grow in mass. 



Table 2. Properties of helium normal (Yq = 0.24) stars. Masses are in solar 
units, lifetimes in Myr. 



M 


1"MS 


TRGB 




tagb 


M c ,,gb 


Mc.cAGB 


0.60 


37601 


2751.1 


185.0 


7.70 


0.47 


0.38 


0.71 


18999 


2264.3 


109.4 


5.00 


0.49 


0.42 


0.79 


11934 


1976.6 


107.0 


4.30 


0.49 


0.43 


0.89 


8414 


722.7 


107.2 


4.55 


0.49 


0.43 


1.00 


5042 


1019.7 


106.8 


4.20 


0.49 


0.44 


1.10 


3610 


782.5 


104.5 


3.56 


0.49 


0.45 


1.20 


2920 


399.9 


103.8 


3.42 


0.48 


0.45 


1.30 


2303 


210.0 


103.3 


3.64 


0.48 


0.46 


1.40 


1831 


122.0 


102.4 


3.10 


0.48 


0.46 


1.50 


1516 


72.8 


105.4 


3.36 


0.47 


0.46 


1.60 


1277 


45.6 


107.5 


2.87 


0.47 


0.47 
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Table 4. As Table|2]for Y = 0.40. 



M 


TMS 


TRGB 


1"HB 


tagb 


M C ,,GB 


M:,eAGB 


0.60 


12336 


1726.1 


168.5 


7.60 


0.46 


0.39 


0.71 


6438 


1009.2 


105.8 


3.43 


0.46 


0.46 


0.79 


4220 


714.8 


100.9 


2.68 


0.46 


0.49 


0.89 


2644 


657.0 


103.0 


2.26 


0.45 


0.51 


1.00 


1938 


294.9 


100.1 


1.99 


0.45 


0.53 


1.10 


1485 


156.8 


98.3 


1.82 


0.45 


0.55 


1.20 


1174 


84.5 


99.4 


1.71 


0.44 


0.56 


1.30 


940 


53.6 


103.2 


1.58 


0.42 


0.57 


1.40 


776 


33.7 


102.5 


0.73 


0.39 


0.58 


1.50 


646 


23.9 


119.6 


0.63 


0.35 


0.44 


1.60 


547 


17.8 


111.6 


1.28 


0.34 


0.61 



Table 3. As Table|2]for Y = 0.32. 



M 


1"MS 


TRGB 


THB 


tagb 


M c ,tGB 


^CcAGB 


0.60 


21986 


2389.3 


184.8 


8.00 


0.47 


0.38 


0.71 


11072 


1745.1 


110.7 


4.20 


0.48 


0.43 


0.79 


7080 


1338.5 


106.7 


3.57 


0.47 


0.45 


0.89 


4711 


871.7 


105.6 


3.11 


0.47 


0.46 


1.00 


3066 


711.4 


103.9 


2.89 


0.47 


0.48 


1.10 


2345 


368.7 


103.4 


3.02 


0.47 


0.48 


1.20 


1810 


212.0 


102.1 


2.56 


0.47 


0.49 


1.30 


1443 


125.4 


102.6 


2.68 


0.46 


0.49 


1.40 


1195 


65.1 


104.8 


2.26 


0.45 


0.51 


1.50 


981 


48.1 


114.9 


2.23 


0.43 


0.51 


1.60 


838 


29.7 


128.2 


0.85 


0.39 


0.51 



Material of the appropriate envelope composition is accreted on 
the star until its mass matches that of the last pre-flash model. The 
helium core is then allowed to grow to the desired mass due to 
hydrogen shell burning until the helium core mass too matches the 
mass in our pre-flash model. This procedure gives us a core helium 
burning star of the correct total mass, core mass and composition, 
but we ignore any mass loss or mixing that may arise as a result 
of the helium flash itself. This is not an unreasonable assumption 
(Harm & Schwarzschild 1966; Dearborn et al. 2006) although some 
authors have found mixing during calculations of the helium flash 
at very low metallicity ( Fujimoto et al. 1990 , Campbell & Lattanzio 
|2008|>. Recent 3D hydrodynamical calculations of the helium flash 
by Mocak et al. ( 2009 ) are consistent with this picture, but suggest 
that the structure of the star on the approach of the helium flash 
may be different than earlier calculations predict. 

We make no attempt to follow thermal pulses during the AGB 
phase and most of our models terminate at the beginning of the 
tp-AGB. 



3 THE STRUCTURE AND EVOLUTION OF HELIUM 
ENHANCED STARS 

Before we consider the evolution of mergers resulting from colli- 
sions involving helium enhanced stars, we review the properties of 
helium-rich stars and compare them to helium normal (Y = 0.24) 
stars. These results are summarised in Tables[2]-[4] 

These tables give the lifetimes on the main sequence (r MS ), the 



Y = 0.24 
Y„ = 0.32 




4.05 4 3.95 3.9 3.85 3.8 3.75 3.7 

logio T eff /K 

Figure 1. Evolution tracks for a 1.2M Q helium-normal star (solid track) and 
a 1.2M (upper dashed track) and a 1.05M o (lower dashed track) Yq = 0.32 
star. 



red giant branch (t rgb ), the horizontal branch (t hb ) and the AGB 
(t agb , not including the tp-AGB). We also give the core mass at the 
tip of the giant branch (M c1 gb) and at the beginning of the AGB 

(M c , eA GB)- 

Perhaps the best-known property of helium-rich stars is that 
they are hotter and brighter than helium-normal stars of the same 
mass (Figures [T] and [2). In this sense helium-rich stars can also 
masquerade as more metal-poor stars, which are also hotter and 
brighter than metal-rich stars. This has been used to explain the 
existence of a blue metal rich sequence in a> Cen ( [Piotto et al.|2005| > 
and more recently to explain the split subgiant branch in NGC 1851 
( |Han et al.|2009) . 

A direct consequence of their higher luminosity is that a 
helium-rich star of a given mass has a considerably shorter main 
sequence lifetime than a helium normal star of the same mass. A 
helium-rich star with the same lifetime as a helium normal star will 
be less massive. The evolution tracks of such "coeval" stars are also 
much closer in terms of colour and luminosity (Figures[TJand[2]l. In 
other words, although it is true that helium-rich stars are brighter 
and bluer than normal stars for the same mass, this effect is much 
less significant when comparing populations of stars with the same 
age. This is the reason that a helium-rich population in a cluster 
shows up as a broadening (or splitting) of the main sequence rather 
than as an extension of the main sequence. 

Helium-rich stars mimic helium normal stars of higher mass in 
other ways. For Y = 0.24, a convective core appears on the main 
sequence above M = 1.1M Q , while for Y = 0.32 an M = 1M Q 
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Y„ = 0.24 
Y n = 0.40 




4.05 4 3.95 3.9 3.85 3.8 3.75 3.7 

logio T e ff/K 

Figure 2. As figure ^ for Yq = 0.4. The lower dashed track is for a 0.9M e 
star. 



solar mass star already has a convective core. For Y = 0.4 the 
convective core already develops for M = 0.9M G . 

The mass loss prescription {T} predicts a higher mass loss rate 
on the RGB for hotter stars, which means that the mass loss rate of 
helium enhanced stars is higher than that of helium normal stars. 
Nevertheless, the integrated mass loss on the RGB is lower for he- 
lium enhanced stars. This is due to the lower luminosity at the tip 
of the RGB, which is in turn due to the helium flash occurring at 
lower core mass. This reduces the integrated mass loss mainly be- 
cause the mass loss rate is highest near the tip of the RGB and also 
because it shortens the RGB lifetime. 

The lower core mass for helium ignition is due to an increased 
temperature at the base of the hydrogen burning shell. The easiest 
way to see why a higher helium abundance in the envelope leads 
to an increase of temperature in the core is by shell source homol- 
ogy ( |Refsdal & Weigert|1970||Kippenhahn & Weigert|l990] >. The 
temperature at the bottom of the hydrogen burning shell scales as 



T'shcll 



M c 



(2) 



where the molecular weight fi is a suitable "typical" value in the 
burning shell. Because the core is degenerate there is a relation be- 
tween M c and R c , say R c oc M~ a with a > (a = 1/3 for cold 
degenerate matter). Using primes to indicate a helium enhanced 
model, we have 



(3) 



so for the same core mass the helium enhanced model will have 
a higher temperature at the bottom of the burning shell (// > ji), 
which implies a higher temperature in the core. An alternative ap- 
proach that leads to the same conclusion is to apply the virial the- 
orem for the core and consider the effect of changing the composi- 
tion of the envelope (c.f. Kippenhahn & Weigert 1990, § 30.5). 

The lower mass loss on the RGB increases the mass on the 
horizontal branch for helium-rich stars of a given initial mass, but 
because of the difference in lifetime a coeval helium-rich horizontal 
branch star will actually be less massive than a helium normal star. 

The core mass at the beginning of the tp-AGB is higher for 
helium enhanced stars of the same mass but also for helium en- 
hanced stars of the same lifetime. Because of this we expect the 
white dwarf remnants of helium enhanced stars to be more massive 
than the remnants of helium normal stars. This is an interesting re- 



sult but the implications of this for the white dwarf population in 
globular clusters are beyond the scope of this work. 



4 STRUCTURE OF THE COLLISION PRODUCTS 

4.1 Lifetime and luminosity and colour functions 

The structure and evolution of the collision products agrees very 
well with previous findings (Sills et aT]|1997| |Glebbeek & Pols] 
2008 ), bearing in mind the differences between helium enhanced 
and helium normal stars discussed in the previous section. 

As can be expected from the "entropy sorting" principle, the 
most helium-rich object will settle at the core of the collision prod- 
uct. Figure [3] shows the difference in the evolution tracks for an 
M = 0.6 + 0.4M Q collision product with different helium abun- 
dances in the primaries: as long as at least one of the stars has 
Y = 0.24, the evolution tracks look qualitatively the same, but 
brighter and bluer with increasing helium abundance in the inte- 
rior. 

The evolution track of a collision product resulting from Yq = 
0.24 and Y = 0.32 progenitors is intermediate between the evo- 
lution tracks of stars with the same mass but an initially uniform 
helium composition of Yq = 0.24 or Yq = 0.32. For parent stars of 
equal mass the evolution track is similar to that of a Y = 0.28 star 
of the same mass although it is not as blue. The evolution of the 
core, and therefore the lifetime, more closely resembles that of a 
Y = 0.32 star than that of a Y = 0.28 star. 

In previous work (Sills et al.||1997| |Glebbeek &Toisl[2008) 
we showed that the core of the collision product is formed by the 
secondary, unless the primary has evolved sufficiently to reduce 
its central entropy below that of the secondary. This is slightly 
more complicated in the case of merging stars with different he- 
lium abundances. An unevolved star with a higher helium abun- 
dance will have a lower central entropy than an unevolved star of 
the same mass with a normal helium abundance. It will therefore 
tend to settle in the core of the collision product, which will resem- 
ble a more evolved star, or the collision product of more evolved 
parents. For instance, the result of an Mi = 0.6M Q , Yq = 0.24 and 
an M 2 = 0.6M o , Y a = 0.32 collision at 6000 Myr resembles the 
collision product of Mi = 0.4M o , Y = 0.24 and M x = 0.8M o , Y Q = 
0.24 at 8 000 Myr. 

In general, it is hard to predict which star will form the core 
of the collision product without looking at the structure of the two 
stars. It is still true that if the primary is sufficiently evolved, it 
will provide the core of the collision product even if the secondary 
started out with a higher helium content. However, it will need to 
be more evolved than in the situation where the secondary has the 
same helium content initially. 

In either case, the lifetime of the collision product is reduced 
because its central helium abundance is higher than the situation 
where both stars had a normal helium abundance. This effect is ex- 
acerbated by the increase in the stars' luminosity due to the higher 
helium abundance in the envelope. 



4.2 Comparison with observed blue straggler populations 

We compared the results of our models with the blue straggler pop- 
ulations in four galactic globular clusters: NGC 2808, NGC 1851, 
5634 and NGC 6093 (M80), see Table[5] The latter cluster was also 
studied by Ferraro et al. (2003). Observations were taken from the 
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Figure 3. Evolution tracks for collisions between an M = 0.6M o and an 
M = 0.4M e star at 6000 Myr. The initial contraction phase is not shown for 
clarity. 



Table 5. Photometric parameters of our comparison clusters, taken from 
|Harris|(l996) . 



ID 


Name 


E(B - V) 


(m-M) 


[Fe/H] 


NGC 1851 




0.02 


15.47 


-1.22 


NGC 2808 




0.23 


15.56 


-1.15 


NGC 5634 




0.05 


17.16 


-1.88 


NGC 6093 


M80 


0.18 


15.56 


-1.75 



HST WFPC2 database by |Piotto etaLlpOO^ and blue stragglers 
were selected using the procedure of |Leigh et al.|p007| lp] 

These clusters were chosen because they all show different in- 
dications for multiple populations. In the case of NGC 2808, there 
is an observed splitting of the main sequence in at least three se- 
quences ( |Piotto et al.|20 07 1, which is attributed to the cluster popu- 



Blue straggler numbers for the cores are given in |Leighetal.|j2008) . We 
use the same procedure to find blue stragglers in the outer regions of the 
clusters. We find 33, 50, 11 and 17 blue stragglers for NGC 1851, 2808, 
5634 and 6093 respectively. 




Figure 4. Colour-magnitide diagram for NGC 2808, showing the blue strag- 
glers and a 12 Gyr isochrone for our adopted E(B - V) and m-M (solid 
line) as well as the E(B - v) = 0.23 (dashed line). 



lation consisting of distinct sub-populations with a different helium 
abundance. It has become common to refer to these as "first gener- 
ation" and "second generation" (or "third generation"), where the 
"first generation" refers to stars with a normal (low) helium abun- 
dance and "second generation" refers to the helium enriched popu- 
lation that (supposedly) formed later than the "first generation" out 
of polluted material. 

For NGC 1851, the sub-giant branch is known to be split in 
at least two distinct populations (Sav iane et al.|19 98 ; Milon e~et al.| 
|2008| >, which can be explained if there are two populations with 
different total CNO abundances (Cassisi et al. 2008, Ventura et al. 



2009) but with the same helium abundance. Recently, Han et al. 
(2009) have presented evidence for a split giant branch in U — I 
colours. Their interpretation of this is that there is a Helium en- 
hanced population in NGC 1851 (with Y = 0.28) that also has a 
slightly higher metallicity and they show that this would only show 
up as a split horizontal branch in V - I colours. 

Apart from the usual globular cluster abundance anomalies 
there is no photometric indication for a "second generation" in 
NGC 6093 or NGC 5634. 

The scenario where some blue stragglers are formed through 
collisions between stars of different helium content should there- 
fore apply only to the case of NGC 2808. 

Comparing the photometry from Piotto et al.| (2002[ l to that 
of |Walker| j 1 999 ), there appears to be a systematic colour shift 
A(fi - V) « 0.1 to the blue in the |Piotto efaT] < f2002] ) data (be- 
fore applying reddening corrections). The origin of this offset is 
unclear, but unless we correct for it we are unable to properly fit 
an isochrone to the data. In practice, we adopted the Piotto et al. 
(2002 ) data as it is available and determined an effective value of 
E(B - V) = 0.07 for our isochrone. We use this value when com- 
paring our models to the Piotto data. 

The reddening for NGC 2808 is somewhat uncertain. Walker 
(1999) finds an overall value E(B - V) = 0.20 ± 0.02, but notes 
that |Schlegel et al.| ( |1998} extinction maps indicate differential red- 
dening across the cluster. Given the offset between the Piot to et al.| 
( 2002 ) and the Walker ( 1 999 ) photometry and the error bars on the 
reddening, our effective value E(B - V) = 0.07 for the Piotto data 
appears to be reasonable. 

In order to compare our models with observations we calcu- 
late the colour and luminosity functions predicted by our models at 
an age of 12Gyr. The mass spacing of our model grid is not very 
fine and this results in very clear discrete jumps and gaps in the the- 
oretical colour and luminosity distributions. To produce smoother 
distributions, we computed evolution sequences for collisions be- 
tween stars with 3 1 distinct masses between 0.4 and 0.8 (inclusive). 
For the masses that fall between those in our grid we interpolate 
between the neighbouring evolution tracks in a similar way to |Pols| 
|etaLl ( fT998) . 

In our model sets collisions only happen at discrete time in- 
tervals. In reality, collisions may happen at any time and for our 
simulated population we should generate collision models at inter- 
mediate times. We account for this by selecting for each evolution 
sequence the portion of the evolution track that falls within an age 
range of 12 Gyr ± 500 Myr. Every point along the evolution track is 
then assigned a weight 



Wi = A/,0(mi)<^(m 2 ) l F(mi,m2,fcoii), 



(4) 



where Af, is the amount of time spent in the vicinity of the current 
point, 0(m) oc rrT" is the initial mass function of Kroupa ( 2001 ) 
and f(mu m 2 , t con ) is the collision probability for star 1 and star 
2 at the time of collision f co n, which we take to be constant for 
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B - V 

Figure 5. Theoretical colour-magnitude diagram for Yq = 0.24 and Z = 
0.001 at 12 Gyr. The solid line is a single star isochrone. 

simplicity. For model sets D and E the relative abundance of the 
different populations is important. We have taken all populations to 
be equally abundant, so that for model set D mixed Y Q = 0.24 and 
Yq = 0.32 collisions are as likely as collisions between stars with 
the same helium content. 

The evolution tracks are then binned in the B— V vs. M v plane 
with the value of each bin being the sum of the weights Wj of evo- 
lution tracks that pass through it. The width of our bins is 0.1 in 
M v and 0.089 inB — V. To convert between theoretical log g, log L 
and T e g to observational My and B - V we made use of the spec- 
tral library by Lej eune et al.| (T997 19981. The resulting colour- 
magnitude diagram for Y = 0.24 is shown in Figure [5] along with 
a Yo = 0.24 isochrone. There are some gaps in the shading in the 
Hertzsprung gap between the main sequence and the giant branch 
where our binning method misses the stars as the evolve quickly 
through this region. Because the time spent in this region of the 
colour-magnitude diagram is small anyway, the weight Wj at these 
points is also small and the presence of these gaps does not affect 
the results of our comparison. 

The collision products detach from the zero-age main se- 
quence around B — V = 0.5 and fill up the blue straggler region 
above and to the blue of the main sequence turn-off. The collision 
product giant branch appears just to the blue of the normal single 
star giant branch and the horizontal branch is slightly bluer and 
brighter than the normal single star horizontal branch. 

Observationally, blue stragglers are selected by drawing a se- 
lection box in the colour magnitude diagram, with cut-offs at the 
low luminosity end to separate the blue stragglers from the main 
sequence and at the high luminosity end to separate the blue strag- 
glers from the horizontal branch. See Leigh et al. (2007} for a de- 
scription of an algorithm to define these selection boxes consis- 
tently between different clusters. In some cases (most notably NGC 
6093 and NGC 1851) the resulting selection box includes one or 
two blue stragglers that are substantially (> O.ldex) fainter than the 
rest. This strongly affects the shape of the luminosity function at 
the low-luminosity end and we have redrawn our selection box to 
reject these stars from both the observations and our models. 

To compare our models to the observations, we have selected 
the blue stragglers from our models using the same selection box 
as was used to select the observed blue stragglers. These selec- 
tion boxes vary slightly from cluster to cluster, so that the range 
of colours and luminosities spanned by the blue stragglers also dif- 
fers slightly from cluster to cluster. 

We compare our Z = 0.001 ([Fe/H] = -1.3) models with 
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Figure 6. Cumulative luminosity functions ip for the observed blue straggler 
populations in the whole cluster (thick red solid line), in the cluster core 
(thick red dashed line) and the luminosity function for the collision models 
(thin blue lines). The solid line is for model set D, the dotted lines are for 
Y = 0.40 (model set C), Y = 0.24 (model set A) and Y = 0.32 (model set 
B) respectively from left to right. 
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Figure 7. Total B-V colour distributions of the observed blue straggler 
populations (histogram) and the colour distributions for the collision models 
(curves). The solid line represents model set D. The dashed, dotted and 
dash-dotted lines represent model sets A - C respectively. 



Table 6. K-S test probabilities comparing the theoretical and observed blue 
straggler luminosity functions for the entire cluster. 
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NGC 185 1 and NGC 2808, which have a slightly higher metallicity 
(Table[5](. Our Z = 0.0003 ([Fe/H] = -1.82) models are compared 
with NGC 5634 and NGC 6093, which straddle this metallicity. 

The (observed and model) blue straggler luminosity functions 
are shown in Figure [6] and the colour distributions are shown in 
Figure|7] Table|6]gives the probability that the observed luminosity 
function is drawn from model sets A, B, C, D or E (defined in Table 
fill, as determined by a standard K-S test. The listed value p is the 
probability using the observed value for M - m. The value p' was 
derived by treating the distance modulus as a free parameter that 
was then derived from fitting to the data. The difference between 
the observed value of M - m (listed in Table [5} and the best fitting 
value for M — m is listed as A(M - m). We will discuss each of the 
clusters in turn. 



4.2.1 NGC 1851 

For NGC 1851 most of our model luminosity functions fit the ob- 
servations about equally well (p ~ 0.45) with the exception of the 
helium enhanced Y = 0.32 models (model set B). If we allow for a 
shift in M - in then the best fits are for model set C (p = 0.90) and 
A (p = 0.75). 

The high luminosity end of the distribution best fits the pure 
Yo = 0.24 models (model set A), with p > 0.999 if we only select 
blue stragglers brighter than My = 2.7. We can match the fainter 
end of the luminosity function, which then matches best with the 
Y = 0.4 collision models (model set C) with p = 0.81. However, 
in that case our models do not match the high luminosity end of the 
luminosity function. The luminosity function for the core is clearly 
brighter than for the cluster as a whole, which could be a signature 
of mass segregation. However, the number of stars is small so it is 
hard to draw any firm conclusions. 

The predicted colour distribution from our models gives a 
tolerable fit to the observations. We predict more stars between 
B-V = 0.2 and 0.3, but the number of observed stars in this colour 
range is small. The predicted location of the peak in the colour dis- 
tribution may be somewhat too blue. A K-S test for the cumulative 
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colour distribution functions is not very conclusive. Model set C 
is too blue on average and has the lowest probability, of 0.52. The 
next lowest probability is 0.97 for model set E. All other model sets 
have probabilities > 0.99. 

Taking into account both colour and luminosity, model set A 
(Y = 0.24) is most consistent with observations for NGC 1851. 



4.2.2 NGC 2808 

For NGC 2808 the best agreement (p = 0.62) is obtained for a 
mixed population of collision products involving Y = 0.24 and 
Yo = 0.32 stars (model set D), although good agreement with the 
pure Y = 0.24 and a mixture of all collision products can also be 
obtained by allowing for a shift in the distance modulus. The com- 
bination of Yo = 0.24 and Yq = 0.32 marginally remains in best 
agreement with the observations, however. Interestingly enough, 
the high-luminosity end of the luminosity function is in better 
agreement with the pure Yq = 0.32 collisions, which is in line with 
expectations. 

The core luminosity function seems to fit slightly better with a 
pure Y = 0.24 component for luminosities that are one magnitude 
brighter than the turnoff, which is compatible with mass segrega- 
tion since stars of Y = 0.24 are expected to be more massive than 
stars with Y = 0.32 on average. 

The range of colours spanned by our collision models agrees 
well with the observations, and both place the peak between A(B - 
V) = -0.2 and -0.1. The observations may fall off more quickly on 
the blue side, but the difference is within the expected error based 
on the number of observed stars in the colour bins. On the other 
hand, the observations clearly show an excess of blue stragglers to 
the red of A(B - V) = -0.1. In our models the stars in this colour 
range are post-main sequence objects that are in the Hertzsprung 
gap. Similar discrepancies between models and observations have 
been noted before (e.g. Sills et al. 20001. Extra mixing, for instance 
due to rapid rotation, offers one possible way to extend the lifetime 
of stars in this region. Convective overshooting, as noted before, is 
ineffective in removing this discrepancy. Another possibility is that 
at least some of the stars in this region are unresolved binaries. 

A K-S test of the cumulative colour distribution gives about 
equal probability of 0.70 to both model sets B and D, followed by 
model set A with a probability of 0.62. The pure Y = 0.4 models 
(model set C) has a probability < 0.001. Combining both the colour 
and luminosity information, our model set D best describes the blue 
stragglers in NGC 2808. 



4.2.3 NGC 6093 

None of our luminosity functions fit particularly well. The best fit- 
ting model is model set A, although model set D is not much worse. 
By varying the distance modulus, model set A can be made to fit 
a bit better, while the fit for model set D does not improve signifi- 
cantly. 

The colour distributions for the different model sets are all 
very similar. The observations show a peak that is 0.15 magnitudes 
bluer than the tumoff. As with NGC 2808, there are more observed 
blue stragglers redward of A(6 - V) = -0. 1 than predicted by our 
models. By contrast, no blue stragglers are observed blueward of 
A(S- V) = -0.3. The Y = 0.24 models fit slightly better (p = 0.39) 
than the mixed Y = 0.24; 0.32 models (p = 0.25) but the picture is 
not very clear. 



4.2.4 NGC 5634 

Again, none of our models fit very well. Model set A (Y = 0.24) 
gives the best agreement with the observations, although model set 
D is also not a bad fit if we vary the distance modulus. The theo- 
retical luminosity function falls off perhaps a little too quickly for 
higher luminosities. The colour distributions are not very different 
for the three model sets and all match about equally well. The ob- 
servations show a peak 0.1 dex bluer than the turnoff that is not 
present in the models. 



4.3 Late evolutionary phases 

The late (post-main sequence) evolution of collision products is in- 
teresting for two reasons. First of all it allows us to test our under- 
standing of the subsequent evolution of collision products by com- 
paring the observed distributions and properties of evolved blue 
stragglers (post-blue stragglers) with the models. Second of all, it 
may allow us to probe the dynamical history of the cluster over a 
longer time interval than can be accomplished with the blue strag- 
glers alone. 

The most promising post-main sequence evolutionary phase 
to identify stellar collision products is during core helium burn- 
ing. Observationally, the collision products are then expected to lie 
above the horizontal branch ( Renzini & Fusi Pecci 1 988 ; Fusi Pecci 
|et al.|1992||Ferraro et al.|1999||SiHs et al.|2009| >. 

We can derive a selection box for evolved collision products 
by using the evolution models. The first step is to determine the 
minimum and maximum values of My and B - V during core he- 
lium burning for all of our collision models. These define a selec- 
tion box in the colour-magnitude diagram that contains the core 
helium burning phase of all our collision products. However, for 
some models this selection box will now encompass more than just 
the core helium burning phase. Due to the difference in helium con- 
tent, the colour of the giant branch is different for different models 
in our set. We therefore restrict our selection box so that it does not 
overlap with the red giant branch in any of our models. We also 
impose a cut to remove the cluster horizontal branch. This way, 
we aim to select only stars that are unambiguously evolved blue 
stragglers. Because we have restricted the selection box, we will 
not capture all of the core helium burning phase for all collision 
products, and we may miss ~ 20% of the core He burning colli- 
sion products by using this narrower selection box. In the log 10 T s g 
I log 10 L/L e plane our selection box is very similar to that of Sills 
[etaLlp509) . 

We can define a selection box for the AGB in the same way 
as for the core helium burning phase, but this selection box turns 
out to be very narrow. Since the expected ratio of blue stragglers to 
AGB post-blue stragglers is also very high, ~ 1000, it is virtually 
impossible to compare with observations. Therefore, we only com- 
pare the ratio of blue stragglers to core helium burning post-blue 
stragglers. 

The predicted ratio of blue stragglers to evolved blue strag- 
glers is typically 20 - 50. Observationally, the ratios are about 10, 
but there are only a few stars in the evolved blue straggler selection 
box and numbers are very sensitive to whether any particular star 
is classified as an evolved blue straggler or not. 

A reasonable agreement between models and observations 
was reached by |Sills et al.| (|2009 ) based on the average of the hor- 
izontal branch and main sequence lifetimes of their collision prod- 
ucts, as long as they selected only the brighter (more massive) col- 
lision products. Our approach differs from theirs in two respects: 
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first of all we calculate a population of collision products where 
each collision product is assigned a weight that depends on the IMF 
probability of the parent stars as opposed to taking a straight aver- 
age. Second of all we calculate our population ratios, for both the 
models and the observations, based on selection boxes in the colour 
magnitude diagram, not on the relative duration of the evolution 
phases themselves. 

The first of these has the effect of giving more importance to 
the lower mass collision products, which will tend to increase the 
predicted ratios. The effect of using selection boxes that will only 
capture part of the evolution rather than comparing the lifetimes di- 
rectly will similarly tend to increase the blue straggler to post-blue 
straggler ratio. Neither of these effects should affect the compar- 
ison between our models and the observations, however, because 
blue stragglers and evolved blue stragglers are selected consistently 
between the observations and in the models - if the models reflect 
the observed population, then these population ratios should come 
out the same as long as the selection criteria are the same. How- 
ever, the small number of observed stars that actually fall within 
our selection box make it hard to draw any firm conclusions from 
this comparison. 



5 SUMMARY AND DISCUSSION 

We present new evolution calculations for stellar collision prod- 
ucts, where we have allowed the progenitor stars to have a different 
helium abundance. 

Our collision models do a reasonable job of reproducing the 
observed blue straggler luminosity function and colour distribution. 
We predict, perhaps, too few blue stragglers at the red edge of the 
observed distribution. For NGC 1851, the best agreement is ob- 
tained for a single population of helium normal stars (Y = 0.24) 
while for NGC 2808 the best fit is obtained with a population of 
mixed Yq = 0.24 and Yq = 0.32 stars. These results are what we 
would expect in light of observations of multiple populations in 
these clusters. A lower metallicity set of models agrees with NGC 
6093 and NGC 5634 in the sense that the best fitting models are 
those without a helium enhanced population. The observed popu- 
lation ratios of blue stragglers to evolved blue stragglers is larger 
than is seen in the observations, but the number of observed evolved 
blue straggler candidates is small. 

Recently, |Han et aT] f 20091 reported the presence of two dis- 
tinct populations in NGC 1851, one helium normal and one helium 
enhanced (Y = 0.28). However, they also report that the helium- 
rich population has a higher metallicity, which offsets the colour 
changes due to helium enhancement, at least in V — I colours. Since 
all our models have the same metallicity, there is no corresponding 
compensation for the shift in colour due to helium enhancement. 
To compare our models more directly we would need to allow for 
a similar shift in metallicity. 

Although the agreement between our models and the obser- 
vations is reasonable, the agreement is not perfect and the models 
presented here are not complete. There are two obvious improve- 
ments that can be made. The first of these has to do with the blue 
straggler formation rate, which we have effectively taken to be con- 
stant and independent of the helium abundances of the two col- 
liding stars. By allowing the collision rate ¥ in j4]l to vary with 
'con we have more freedom in shaping the colour and luminosity 
functions - in fact, if we want to use blue stragglers (and possi- 
bly evolved blue stragglers) to learn about the dynamical history of 
their host cluster, then it is essential that we allow this factor to vary. 



Simulations of clusters hosting multiple populations, like those of 
|D'Ercole etaT] ( 2008 1 may serve as a guide for how this should 
be done. Ideally, the collision rate should be directly determined 
from dynamical simulations of cluster evolution. A software envi- 
ronment that might be especially suited for this task is the MUSE 
software package, which is in active development (Portegie s Zwart| 
|et al.|2009| >. 

The second improvement that can (and should) be made is that 
we should consider binaries in addition to single star models. For 
the blue stragglers, this is an obvious addition because binary mass 
transfer is an alternative scenario for blue straggler production (e.g. 
Chen & Han 2009 1. However, binaries are also important in another 
respect: observationally, it is not possible to distinguish the light 
of the two stars in a binary system, which may place the system 
in an unusual point in the colour magnitude diagram that is hard 
to reproduce with single star evolution tracks. This becomes espe- 
cially important when we look for evolved blue stragglers because 
a blend of a horizontal branch star and a red giant may appear in the 
same region of the colour-magnitude diagram as the evolved blue 
stragglers. Such unresolved binaries can be recognised by multi- 
wavelength photometry because the colour of the unresolved bi- 
nary will not change in the same way as that of a single star, and 
the binary will move to a different region of the colour-magnitude 
diagram: a normal star will stay close to other stars of a similar 
spectral type, but an unresolved binary will move closer to the po- 
sition of one of its unresolved components. 

Evolved blue stragglers offer an interesting possibility to test 
our understanding of blue straggler evolution, but because the num- 
ber of observed post-blue straggler candidates is small it is espe- 
cially important to understand how this region of the colour magni- 
tude diagram may be influenced by the presence of binaries. Both 
normal binary interaction and collisions can increase the actual 
number of stars in the evolved blue straggler region without in- 
creasing the number of blue stragglers. Unstable binary mass trans- 
fer from a red giant to an unevolved main sequence star companion 
can lead to a spiral in and merger of the two stars if the envelope is 
not ejected. A collision with a red giant will have a similar result. 
The red giant core is likely to remain intact during the merger so the 
merger product will still have a degenerate helium core and evolve 
like a more massive red giant. Such a merger would show up as an 
"evolved blue straggler" despite never having been a blue straggler 
itself. 

Our population models are consistent with observations of 
multiple populations in the sense that helium enhanced model sets 
fit best with clusters (in particular, NGC 2808) where helium en- 
hancement has been inferred from the observations. We hope that 
with the inclusion of a population of binaries, population models 
such as we have presented in this paper could be used not only to 
test our understanding of cluster dynamics but also to get a better 
handle on the nature of the multiple populations that are now ob- 
served in star clusters. Future cluster simulations that include an 
accurate treatment of the evolution of stellar collision products as 
well as multiple populations will be an important diagnostic tool 
and we plan to make our models available for such a study. 
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